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In the continuum, a topological obstruction to the vanishing of the non-abehan anomaly in 2n 
dimensions is given by the index of a certain Dirac operator in 2n+2 dimensions, or equivalently, 
the index of a 2-parameter family of Dirac operators in 2n dimensions. In this paper an analogous 
result is derived for chiral fermions on the lattice in the Overlap formulation. This involves deriving 
an Index Theorem for a family of lattice Dirac operators satisfying the Ginsparg-Wilson relation. 
The index density is proportional to Liischer's topological field in 2n+2 dimensions. 

o" 

o 

fSJ ' The Atiyah-Singer Index Theorem, both for single operators [Q and families of operators [||, has been of major 

importance in the modern development of Quantum Field Theory. For example, the Index Theorem for single operators 
gives topological insight into the non- vanishing of the axial anomaly H , provides the basis for a resolution of the 
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U(l) problem W and determines the dimension of the instanton moduli space (used in semiclassical investigations of 

^\ Yang-Mills theory; see, e.g., P), while the Families Index Theorem reveals topological obstructions to the vanishing 

CN . of gauge anomalies, thereby providing constraints on allowable theories ||-||]. 

When attempting to get a well-defined non-perturbative formulation of QFT's (in particular, gauge theories) one 

~^ , successful and well-established approach has been to formulate the theory on a spacetime lattice M. Therefore, 

vQ ' carrying over the Index Theorem, both for single operators and families, to the lattice is an interesting and important 

fT^ , problem. For a long time this did not seem possible though, due to the fermion doubling problem and the resulting 

^^ ' need for acceptable lattice Dirac operators to break chiral symmetry pO| . In the traditional formulations, at best only 

^^ , a remnant of the Index Theorem is retained on the lattice pj. However, the situation has changed quite dramatically 

^~j ' in recent years with the advent of the Overlap formulation ]12[| and the discovery p^p^ of acceptable lattice Dirac 

^y. , operators satisfying the Ginsparg-Wilson relation ||l^ 

"t^ I D^5 + 75Z) = aD^^D (a=lattice spacing) (1) 

JL . Such operators have exactly chiral zero-modes (since Dtp = ^ D{'-fc,tp) = [aD^^D — 'yc,D)i/j — 0), which allows to 

qj ■ define index £) = Tr(75|kcrD) @. There is a "Lattice Index Theorem" ||^,^,[8) 

^ ■ 

'^[ indexZ? = -^Tr(75i?)=a4^(?(a;). (2) 



where 



q{x) = --tr(75D(x,x)) (3) 



is the index density. For SU{N) gauge fields on the Euclidean 2n-dimensional torus, index!? and q{x) reduce to the 
continuum index and density in the classical continuum limit |19| , at least when D is the overlap Dirac operator [[l3| 



(Earlier results in this direction were obtained in |15| , |l8 20 1. When D is the overlap Dirac operator the right-hand 



side of (g) has a spectral flow interpretation which had previously been used as a definition of lattice topological 
charge in |12(| .) Furthermore, although it is not invariant under the usual chiral transformations, the fermion action 
S = a'^'^^ip{x)Dtp{x) exhibits an exact lattice-deformed version of chiral symmetry |Tq] (which was implicit in the 
overlap formalism): SS — for 6ip = "/^ip , Stp = ■075 where 

75=75(1-0^) (4) 

An easy consequence of (j^ is 75 = 1. Furthermore, after supplementing (yj) with the 75— hermiticity condition 

D* = 75^75 (5) 

we have 75* =75. Thus 75 can be viewed as a lattice-deformed chirality matrix. The axial anomaly for the lattice- 
deformed chiral symmetry transformation above can be determined from the corresponding change in the fermion 



measure to be A{x) = —iatr{'-f5D{x, x)) — 2iq{x) |16). This is completely analogous to the relation between axial 
anomaly and index density in the continuum l3l . 

Having seen that there is an exact Lattice Index Theorem for lattice Dirac operators satisfying the GW relation, 
and that the index and its density are related to the axial anomaly in precisely the same way as in the continuum, it 
is natural to ask if there is also a Lattice Index Theorem for families of such operators such that the families index 
is related to gauge anomalies (or more precisely, to obstructions to the vanishing of these anomalies) in the same way 
as in the continuum. In this paper we show that this is indeed the case: We derive an Index Theorem ((O) below) 
for a family of such lattice Dirac operators, parameterised by a 2-sphere in the orbit space of SU{N) lattice gauge 
fields on the 2n-dimensional torus T^". This is the prototype for a more general Lattice Families Index Theorem 
which is currently under development [ pi[ . We find that this index is related to an obstruction to gauge- invariance in 
precisely the same way as in the continuum setting, where it was previously studied by Alvarez-Gaume and Ginsparg 
0. Furthermore, the index density is found to be proportional to Liischer's topological field q{x, 1/1,1/2) in 2n-|-2 



dimensions |23| (given by (|26| ) below) . This provides a natural origin for q{x,yi,y2) in the lattice theory. (It was 
introduced in an ad hoc manner in |23|] .) This is of interest and potential use in connection with Liischer's approach 
towards achieving gauge-invariance in nonabelian lattice chiral gauge theory: a local gauge anomaly-free formulation 
exists if and only if the local cohomology class represented by q{x, y2, 3/2) is trivial Q. 

In the continuum, the modulus of the chiral determinant det{i/j^) (suitably regularised as in p|) is gauge-invariant, 
but anomalies may arise in the phase. Consider a family (pg of gauge transformations parameterised by G S^. Each 
(f>e is a map T^" -^ SU{N) (we assume for simplicity that the fermion is in the fundamental representation). The 
action of (pg on A determines a circle-family {A^}g^gi in the space of continuum SU{N) gauge fields on T^". We 
restrict A to be in the topologically trivial sector (otherwise the chiral determinant vanishes). Then, generically, 
det(/J^) 7^ and we have a map 

S^ ^S^ CC , 0<^ det{i/^f)/det{i/^^) (6) 

(where 5^ C C denotes the unit circle in C). The winding number Wc of this map is an obstruction to gauge- 
invariance of the chiral determinant (since if the determinant is gauge-invariant then it is constant around the family 
{A^}g^gi and the winding number is zero). In Q the winding number Wc was shown to equal the index (which we 
define here to be the topological charge of the index bundle) of a 2-parameter family of Dirac operators /D^^'*\ or 
equivalently, the index of a Dirac operator Vc in 2n+2 dimensions, given as follows. (The subscript c here and in the 
following refers to "continuum".) The family A^ is extended to a family A^^'*) = tA^ with {6,t) € B^, the unit disc. 
^|^(s,t)| corresponds to a 2-sphere in the orbit space of gauge fields since the A^^'^^'s are all gauge-equivalent.) This 
determines the family of Dirac operators ^^^'*) =/)'^ ' = 7^(9^ 4- A}^ ). The Dirac operator Vc acting on spinor 
fields on B^ x T^" is now given by 

Vc^V^iidc^+Aa) a = l,...,2n + 2 (7) 

where the gauge field A on B'^ x T^" is given by -4^(^, t, x) — AJj, '*' (x) , Aa = for a = 2n + f , 2n + 2, and the Dirac 
matrices in 2n-|-2 dimensions are chosen as F'' — ai^ 7'' , F^"+^ = cti <8i 75 , F^"+^ = 0-2 I where ctj (j — I, 2, 3) are 
the Pauli matrices and 75 = ^"7^7^ • • • 7^". The 7^'s and ctj's are taken to be hermitian so^ is anti-hermitian and Vc 
is hermitian. In (R) the derivatives for a = 2n -I- 1, 2ri -I- 2 are d2n+i — -^ , d2n+2 — -^ where (yi, j/2) is a cartesian 
coordinate system on B^ . The polar coordinates 9 and t are henceforth viewed as functions of (j/i, 2/2)- Let B^ denote 
another copy of the unit disc, then B^ x T^" and B^ x T^" can be glued together along their common boundary 
S^ x T^" to get the closed manifold S"^ x T^". An SU{N) vectorbundle over this manifold is defined by taking the 
transition function on the common boundary S^ x T^" to be $(0, x)~^, where $(6*, x) = (j)g{x). The topological charge 
(Pontryargin number) of this bundle is then —deg{$), i.e. minus the degree (generalised winding number) of the map 
$ : S"^ X T^" -^ SU{N). The operator Vc above extends in a natural way to a Dirac operator (also denoted Vc) on 
the spinor fields in this vectorbundle [|6| . Vc anticommutes with the chirality matrix 

F5 = i"+^T^r^ ■ ■ ■ F2"+2 = f73 ® 1 (8) 

and thus has a well-defined index. The main result of H is that the obstruction to gauge-invariance discussed above 
is determined by this index: 

Wc = -index 2?c (9) 

The index of Vr can be calculated from the formula 



indexP, = Tr(r5|kcrpJ = Tr(r5 e""^') Vr > (10) 

This can be calculated in the r — > limit by familiar techniques, leading to ||6| 

indexPc = -deg{<^) (11) 

Thus the obstruction is determined to be Wc = deg{^). 

We now put a hypercubic lattice on T^", with lattice spacing a, and proceed to describe a lattice version of the 
preceding. Chiral gauge theory can be formulated on the lattice in the overlap formalism [ p2[ . This approach can 
be reformulated as a functional integral approach with lattice Dirac operator D satisfying the GW relation (|]) and 
75— hermiticity condition (g) ||2^.|2^. Let C denote the finite-dimensional space of lattice spinor fields on T^". The 
chiral projections P± — ^(1 ± 75) and P± = i(l ± 75) determine decompositions C = C+ S) C- and C = C+ Q> C- 
respectively. The (right-handed) lattice chiral determinant in this setting is 

det{iDl) = {v-,w+{U)) (12) 

where V- and w^ are unit volume elements on C_ and C-^- respectively. These are unique up to phase factors; they can 
be written as w_ = ui A • • • A w^ and i/)+ = wiA- ■ ■ Wd where vi, . . .Vd and lii, . . . ,Wd are orthonormal bases for C_ and 
C+ respectively. w_ and u)+ are the many-body groundstates in the overlap formulation |T3], and correspond to the 
chiral fermion measures in the formulation of ref.'s P2,E3|. Note that 75 = 75(1 — aD^) depends on the lattice gauge 
field U , so the subspace C+ and volume element w+ likewise depend on U . On the other hand, since the usual chiral 
decomposition C = C+®C- does not involve U , neither does u_. We are assuming dimC± = dimC± = d (otherwise the 
chiral determinant vanishes). This is equivalent to assuming index Z)'^ = , i.e. U is in the topologically trivial sector 
[0,p3[. The space of lattice gauge fields will typically contain a subset of measure zero where D^ is not defined. In 
the case of the overlap Dirac operator such fields can be excluded by imposing a condition of the form 1 1 1 — [/(p) 1 1 < e 
on the plaquette products of U pq] . This condition is automatically satisfied close to the classical continuum limit 
since 1 — U(j)) — a'^Ff^^{x) + 0{a^. We will assume that the same is true for the general D that we are considering 
here. 

Let {^ejeesi be a family of SU{N) lattice gauge transformations, then the winding number W of the map 

S'^S' , e^{v^,w+{<f>g-U))/{v^,w+{U)) (13) 

is an obstruction to gauge-invariance of the chiral determinant ( |l2| ) just as in the continuum setting. This was recently 
studied in ||24| where W was shown to reduce to Wc in the classical continuum limit. In the following we will show 
that this obstruction is related to the index of a family of lattice Dirac operators, defined as the index of a Dirac 
operator D in 2n+2 dimensions, in complete analogy with the continuum relation (|9|). 

The action of (pg on U generates a circle-family {(fig ■ C/}egsi in the space U of lattice gauge fields. Choose a 
disc-family B^ — {C/'^'*''}(e 4)532 in U such that f/f^'^) = 0^ • [/. (Such a family might not exist in general due to the 
restrictions on U needed to ensure that D is well-defined. However its existence is guaranteed close to the classical 
continuum limit: we can take the lattice transcript of the continuum family A^^'*^'.) This determines a family of 
lattice Dirac operators D^^'*) = D'^ '' . Setting f ^ = cri 1^)75 (where t's = 75^^'*^ is given by (||) with D = L)(^^*)) and 
r^ = r-^ = (72 (8) 1, we define the Dirac operator in 2n+2 dimensions in the lattice setting to be 

V ^ r'^i{da + Aa,) a = 1,2 (14) 

The derivatives are with respect to the continuous cartesian coordinates {y 1,1/2) on B^ and we have introduced a 
continuum SU(N) gauge field A = Aadi/a on B^ with ^^(yi, j/2, 2;) a function of lattice site x as well as {yi,y2)- "D 
extends in a natural way to an elliptic 1st order differential operator on the vectorfields with values in a vectorbundle 
over the closed manifold S'^ = B^ U51 S^ as follows. The fibre of the vectorbundle is C^ (g) C (i.e. the representation 
space of the Pauli matrices tensored with the finite-dimensional vectorspace of lattice spinor fields on T^n-j g^j^j -j.j^g 
transition function at the common boundary S^ of B^ and B'^ is 1 ® $~^ where ^{0) = (fg. A vectorfield in this 
vectorbundle consists of a function '^{0, t) on B^ together with a function '^{9, s) on B^ , both taking values in C^ (g)C, 
and related at the common boundary S^ by 

$(6), 1) = <^{e)-^ ■ *(6i, 1) EE 1 (g) 0-1 • *(6i, 1) (15) 

V is defined on ^ by (|l|), and is defined to act on $ as tfji{da + ^q) where t\j = ai ® 7^5^ and f ^ = F^. The 
gauge-covariance of D implies that £)(^>i) = D'^"'^ = (t>e ° D^ ° 4'g^ and 75' '^' = (f>g o J5 o (jjg^ . Using these it is 



easily checked that T) respects the relation (|T^), and is therefore a well-defined operator on the vectorfields in the 
above vectorbundle over 5^, provided the gauge field A — Aadya on B^ is related to the field A on B^ at the common 
boundary S^ by 

A{0, 1, x) = ^g{x)-^A{0, 1, x)(f>e{x) + (l)g{x)-^dg<j>e(x) (16) 

(E.g. we can take A = Q and A{9, s, x) — s(f>g{x)^'^dg<j>e{x) in terms of polar coordinates {9, s) on B^.) 

The space of vectorfields in the above vectorbundle over S*^ is denoted by V in the following. The chirality operator 
Fs = (73 (g) 1 determines a chiral decomposition V = V+ © V_ . The ellipticity of V follows easily from the facts that 
f 1 ® 75 anticommutes with (72 (8) 1 and (ci (8) 75)^ = ((72 1)^ = 1 (g) 1. Also, I? is formally self-adjoint with respect to 
the natural inner product in V since cti (8) 75 and (72 1 are self-adjoint on C^ ®C. T) anticommutes with Fs = 0^3 (8) 1 

and therefore has a chiral decomposition I ^ „" ] and index T) — dim ker !?_)_ — dim ker I?_ . The following formula 

for the index is derived below: 

Theorem. 

indexP = --^ { f TT(P+dP+dP+) + \ f Tiic^-^de^ei^'^)] (17) 

Here P+ is to be viewed as a function on the space U of lattice gauge fields whose values are operators on C (i.e. 
finite-dimensional matrices), and d is the exterior derivative on hi. Thus the first integrand is a 2-form on hi and can 
be integrated over the disc B^ in U to get a C-numbcr. The second integrand is a 1-form on the boundary 5^ of S^, 
with 75 — 75 '■°'"'^' constant. 

By Eq.(3.11) of EJ] the obstruction (winding number) W associated with the map (O) equals the right-hand side 
of (ll^) without the minus sign. It follows that 

ly^-indexP. (18) 

This is the promised lattice analogue of (ph. Since W reduces to Wc = deg{^) in the classical continuum limit [E4|, 
it follows that index© reduces to index2?c = ^deg{^) in this limit. The 2-form in the first term in the right-hand 



side of (O) has appeared previously in the overlap formalism in |26|, where it was interpreted as a form of Berry's 
curvature. (The Berry phase is associated with the state w+{U) in (p^).) It is interesting to note that a version of 
this 2-form also arises in the context of the quantised Hall effect p^ . The second term in (ll^) arises in p3| as the 
integral of the covariant gauge anomaly (and vanishes in the special case where U = IJ. 

The index formula ( [l7|) , together with ( [25| ) for the index density, and the relation ([l8| ) are the main results of this 
paper. The proof is as follows. Analogously to (|lO|) we have 



index X> = Tr(F5 e^'"^ ) Vr > . (19) 

This can be evaluated in the t ^ limit by the same familiar techniques used to evaluate (|l^ in : It is seen to be 
the sum of a contribution from the B"^ part of S*^ = B^ LI B^, given by 

/ d^ya^y2qr,ix,yi,y2) (20) 

where 

/OC t2 l 

-^ tr(F5 e-''=-^(e--^')e^'=-^)(x, x) , (21) 

and an analogous contribution from the B^ part. In ( [2l] ) the trace is over spinor and flavour indices; 0{x,y) denotes 
the kernel function of an operator O on scalar lattice fields. Setting Vq = da + Aa we find from ( [l4|) that 

_ p2 ^ pay^p/3y^ ^ ^,^/3y^y^ ^ ^a^^ ^ ^^l^^ ^ ^ ^ ^^ ^ ^ ^a375-Fl2 + (75V175 - ^a3V275)Vl 

= (9^ + (2Ai + 75V175 - ifT3V275)ai + 2A2d2 + V„A„ + ^^375^^12 + (75V175 - icr3V27'5)^l (22) 

where Vq75 = dajs + [Aa, 75] (as in ||23| ); for notational simplicity we have omitted the (E) symbol. After substituting 
this in ( pol) and making a change of variables kj -^ T~^^^kj we find 



1 f°° d^k 
gi?(a;,yi,y2) = lim - / tr(g3 exp{-fc^ + ^/t{2Ai + 75V175 - ia'iV2ib)iki + ^{2A2)ik2 

+ t(V„A„ + 75V175 + icJz{l\Fi2 - V275^i))})(2:,2:) (23) 



This can be calculated by expanding the integrand in powers of y^r. Terms with odd powers of CT3 give vanishing 
contribution, as do terms with fc-dependence of the form e~^ fcifef where either p or g is odd. The remaining terms 
are 

rie-'^'tr((2AiV275 +7"5Vi75V275)fc? +75^^12 - ^2l5Ai){x,x) + 0(t^'^) (24) 

where we have used the fact that 75Vq7'5 = — Vc(7'575 (an easy consequence of 75 — 1). Evaluating the integral in 
( P3| ) with this integrand, we find that the contributions from the 2AiV275fcJ and — V275A1 terms in tr(- • •) in ( p^ ) 
cancel, resulting in 

qvix,yi,y2) = 7p-( 7tr(7'5Vi7'5V27'5)(a;,x) + -tr(7'5Fi2)(a;, a;) ) (25) 

Zm V 4 2 / 

Modulo the numerical factor — l/27r, this coincides with Liischer's topological field P3| 

q{x,yuy2) - -^tr(i75[VlP_, V2F^] + i[ViF_, V2P_]75 + iPi27~5)(a:, a;) . (26) 

(Note that 75[ViP_, V2P-] = [ViP_, V2P-]7'5 = 575ViX5V27'5 (an easy consequence of 7''5Vq7'5 = "'^aibjb) and 
that tr(75(a;, a;)Fi2(a;)) = tr(Pi2(x)75(x, a;)).) Summing (|25| ) over the lattice sites gives (cf. Appendix B of P3[ ) 

a'E9^(^'2^i'2^2) = ^Tr(i (75^1755275) - ^^1(^275) + ^^2(^175)) (27) 

The contribution to (EOh from the first term in (B^) gives the first term in the index formula (f^- The contribution 
to ( po| ) from the remaining terms in ( p7| ) reduces to —-^Jgi Tr(A(6', 1)75^ ' ') in polar coordinates. The analogous 
contribution to index I? from the B^ part is only +i J„^ Tt{A{6, 1)75 ) (since 75 is constant). Adding these and 
using ( p^ we get the second term in ( 1^ . 
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